Re-entrant Behavior and Gigantic Response in Disordered Spin-Peierls System 
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Effects of disorder and external field on the competing spin-Peierls and antiferromagnetic states 
are studied theoretically in terms of the numerical transfer matrix method applied to a quasi one- 
dimensional spin 1/2 Heisenberg model coupled to the lattice degree of freedom. We show that, 
at temperatures above the impurity-induced antiferromagnetic phase, inhomogeneous spin-Peierls 
lattice distortions remain to exist showing a re-entrant behavior. This feature can be drastically 
altered by very weak perturbations, e.g., the staggered magnetic field or the change in interchain ex- 
change coupling J^ , leading to a huge response, which is analogous to the colossal magnetoresistance 
phenomenon in perovskite manganese oxides. 

PACS numbers; 75.10.-b, 75.47.Gk, 71.30.-fh, 75.lQ.Pq 

In strongly correlated electronic systems, various kinds 
of long ranged ordered (LRO) phases emerge, frequently 
those energies being very close to each other, i.e., with 
multicriticality. In such cases, the system becomes sensi- 
tive to perturbations such as applied pressure, magnetic 
field, electric field, carrier/impurity doping, etc., leading 
to intriguing new phenomena. 

The disordered spin-Peierls (SP) system realized in 
doped CuGeOa [1] is a typical example showing such sub- 
tle multicriticality of different phases. Its phase diagram 
on the plane of temperature, T, and impurity concentra- 
tion, X, is experimentally explored in detail (Fig. 1(a)), 
which is almost universal regardless of the dopant sub- 
stituting for either Cu or Ge sites [2]. The nonmagnetic 
SP state realized in undoped CuGeOa is fragilely destabi- 
lized by impurities of a few percent, as shown in Fig. 1(a). 
The impurities induce local spin moments where the in- 
terchain exchange coupling results in the antiferromag- 
netic (AF) state, which, at the small- a; and low-T region, 
co-exists with SP lattice distortion in an inhomogeneous 
way (SP-f AF phase in Fig. 1(a)) [3]. On the other hand, 
when X is rather larger, a uniform AF state with homoge- 
neous staggered spin moments becomes stabilized. The 
critical region between the phases is slightly sensitive to 
the dopant, where the system shows either bicritical be- 
havior between the SP and AF states, or tetracritical be- 
havior between all the three phases [2] , though these are 
not thoroughly traced by experiments. Fig. 1(a) shows 
the former case with a bicritical point. Theoretical as- 
pects of this system are also intensively studied analyt- 
ically [4-8] as well as numerically [9], mostly to under- 
stand the novel co-existent state at the small-x region, 
while a few mention the large- a; region [7,8]. 

Another typical example of such interplay between 
multicriticality and disorder in strongly correlated elec- 
tron systems is the colossal magnetoresistance (CMR) 
in perovskite manganese oxides Ri-xA^yhiO^ {R: triva- 
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FIG. 1. Schematic experimental phase diagrams for (a) dis- 
ordered SP system [2] and (b) manganite oxides showing pro- 
nounced CMR [14]. The ordinate is temperature, T, while 
the abscissa represents (a) the degree of disorder realized 
by increasing the impurity concentration and (b) the sub- 
stituent concentration controlling not only the degree of dis- 
order but also the bandwidth (for example Ca concentration 
in Smo.55(Ca,Sr)o.45Mn03 [14]). The grey areas show (a) the 
region of re-entrant SRO of SP state, and (b) the region with 
enhanced CO/00 fiuctuation and CMR effect. 



lent rare earth element. A: divalent alkaline earth el- 
ement) [10]. CMR is seen in compounds where an 
insulator-to-ferromagnetic metal (FM) phase transition 
is observed upon cooling. Small magnetic field can make 
the characteristic temperature such as the maximum in 
resistivity drastically shifted to higher temperature, thus 
resulting in a "colossal" change of the resistivity. To 
understand this effect, the importance of not only the 
multicriticality between the FM phase and the insulat- 
ing charge/orbital ordered (CO/00) state [11,12], but 
also the disorder effect due to quenched chemical dopant 
has recently been pointed out [13,14]. From experi- 
ments carefully controlling the degree of disorder in a 
chemical way [14], the randomness seems to destabilize 
the CO/00 state resulting in the FM ground state [15] 
where the CMR effect becomes prominent, as schemat- 
ically shown in Fig. 1(b). Theoretical works on models 



for manganites with randomness have succeeded to re- 
produce such features [16,17]. 

Here we point out an analogy between the CMR man- 
ganites and the disordered SP system by making a cor- 
respondence of the competing phases as CO/00 ^-> SP 
and FM ^^ AF, respectively. As discussed above, the 
former states are destabilized by impurities, resulting in 
the latter states. As a matter of fact, their phase dia- 
grams on the plane of T and degree of disorder, as seen 
in Fig. 1, are isomorphic to each other. In the disor- 
dered SP system, above the transition temperature for 
the uniform AF state, Taf, short range order (SRO) of 
SP lattice distortion is observed to develop toward Taf 
while it is suppressed below it, resulting in a character- 
istic re-entrant behavior [18]. This is also analogous to 
the CMR compounds where SRO of lattice modulations 
due to CO/00, either static or dynamic, are seen above 
the FM transition temperature [10,14]. The origin of the 
CMR effect is now believed to be due to the sensitivity of 
such re-entrant CO/00 fluctuations [13,14,16,17]. Then, 
similarly, we expect that some "colossal" effect can also 
be seen in disordered SP system, which is investigated in 
this Letter. 

To study the disordered SP system, we consider a quasi 
ID spin 1/2 AF Heisenberg model coupled to the lattice 
distortions [19,20]. The model consists of ID chains, each 
of which being given by the Hamiltonian, 



TiiB ^'^\ji^ + Ui)Si ■ Si+i + — ^u^ 



(1) 



while these are coupled by interchain exchange interac- 
tion between neighboring chains. 






(2) 



where J, J± > 0. Each chain is the usual SP model if Ki 
is uniform, which is unstable toward SP lattice dimcriza- 
tion, Ui = (— l)*Ugp, in order to gain singlet formation 
energy in spite of the loss in elastic energy in the second 
term of Eq. (1) [19]. 

On the other hand, the interchain exchange coupling, 
J±, in Eq. (2) brings about the competing AF state [20], 
to which we apply the mean field approximation as in 
Refs. [8,20]. This is justified for quasi-lD system because 
the correlation length is already large near the transition 
temperature [21]. Then we obtain an effective ID model. 



H' = Hiu zJ^ J2 (i^^)^^ - ^^tfl^) 



(3) 



where z is the number of neighboring chains and the 
mean field approximation for J± is even exact in the limit 
of z — > oo with zJ±^ being finite. Note that in this for- 
malism the mean field {Sf) is independent of chain in the 
original model, and thus consequently ui too [8] . 
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FIG. 2. Theoretical phase diagram for disordered SP sys- 
tem on the plane of T and zj_i_. The phase boundaries for 
the pure case are also written as grey broken lines. The AF 
transition temperature for the pure case, T^p, is identical to 
that in the presence of disorder, Tap for zj± > 0.8. 



Impurities are modeled here by replacing the spin- 
lattice coupling constant, Ki, as Khuik —>■ ^imp ran- 
domly [5,8], and the rcnormalized lattice distortions Ui 
are treated as classical values. We use the numerical 
transfer matrix method [22,23] to calculate the finite-T 
properties of this model, determining Ui and (Sf) for each 
site on the chain self-consistently by minimizing the free 
energy. We use data on open chains with sizes up to 
L = 160 and obtain thermodynamic-limit properties ap- 
plying extrapolation by l/L, and for the random average 
we take 20-40 samples fixing the impurity concentration. 
Note that the method looses accuracy at low-T, at about 
typically T < O.IJ, thus our study is complementary 
to the previous studies on the low-T properties of disor- 
dered SP system [4,8]. In the following we set J = 1 as 
the energy unit. 

In Fig. 2 we show the phase diagram on the plane of T 
and zjj_ for the case of fixed impurity concentration of 
1/8 with K]^uik — 0.8 and Ki„^p — 1.5. In the pure case of 
^"i^buik = 0.8, drawn by grey broken lines in Fig. 2, a bicrit- 
ical behavior is seen at zJi = 0.8 between the SP state 



and the AF state with (Sf ) = (— l)*m^p, rn^p being the 
staggered spin moment, naturally understood from the 
previous studies for the ground state [20]. Here the SP 
transition temperature, Tgp, is constant and the AF tran- 
sition temperature, T^p, is linear as a function of zJ±, 
since we treat Ui as classical values and the interchain 
interaction within mean field approximation. When the 
impurities are introduced, as can be seen in Fig. 2, the AF 
state is extended over a certain range of zJ±. However, 
above this uniform AF phase, we observe the SP state 
and a crossover region in between, where a co-existent 
state of SP and AF with both order parameters spatially 
varying is stabilized, represented as SP-f AF in Fig. 2. 
We shall see this in more detail in the following. 

As an example, in Fig. 3 the evolution of spatial pat- 
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FIG. 3. Spatial pattern of Ui and (Sf) for a certain sample 
of L = 160 with zJi_ = 0.72. Only the middle half of the 
chain (L/4 < i < 3L/4) is shown, to omit open boundary 
effects. The impurity positions with -K'imp = 1-5 (in the bulk, 
^buik = 0.8) are shown by arrows. 



tern of Ui and {Sf) is shown for a certain sample of 
L = 160 for zJ± = 0.72. When T is decreased from 
high temperatures, the SP lattice dimerization is de- 
veloped from slightly below Tgp, but the amplitude is 
modulated such that \ui\ show local minima at the im- 
purity bonds, since Ki^ip > ^buik dislikes lattice dis- 
tortion. When T is lowered further, AF staggered mo- 
ment starts to emerge, locally nucleated by the impu- 
rities out of the SP background, thus co-existing with 
SP in an inhomogeneous way. We note that although 
the amplitudes are modulated, the phases of both SP 
and AF order parameters are unchanged throughout the 
chain [4,8], i.e., there is no solitonic domain wall struc- 
ture as discussed in Refs. [5,6,9]. The co-existence only 
survives at intermediate-T and finally the uniform AF 
state is stabilized at low-T where Ui becomes zero. This 
co-existent state is similar to what is observed in slightly 
doped CuGeOa at low-T [2,4] (see Fig. 1(a)). 

In Fig. 4, the infinite-L extrapolation of the sam- 
ple averaged SP lattice dimerization and AF spin 
moment, -usp = ( x |l]i(~l)*"j| )av and mAF = 
( ■z\J2i(~^y{^i)\ )av, respectively, are plotted as a 
function of T for several values of zJ±. In the region 
of zJ± shown in Fig. 4, first usp develops as T is low- 
ered where the inhomogeneous SP state emerges, while 
toaf = 0, as in the sample in Fig. 3. Eventually ttiaf 
becomes finite and both orders are stabilized as LRO. 
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FIG. 4. Temperature dependence of usp and rriAp for dif- 
ferent values of zj±. Error bars are due to distribution of 
results from samples with different impurity positions. 



There, as ttiaf develops, msp is suppressed, and finally 
the uniform AF state dominates over the SP state. Al- 
though the error bars due to random sample average are 
rather large, we can clearly see the re-entrant behavior in 
Usp, the region with co-existence, and the stability of the 
uniform AF state at low temperatures, to draw the phase 
diagram in Fig. 2. The slopes of the phase boundaries 
of the impurity induced uniform AF state and the co- 
existent state are considerably steep, therefore the states 
can be shifted drastically with small change of the inter- 
chain coupling. 

Next we apply staggered magnetic field Hg by adding 
the term, Hs^^{—iySf, to the effective ID model in 
Eq. (3). We consider this term to control the relative 
energies of SP and AF states, analogous to the uni- 
form magnetic field in CMR manganites. In Fig. 5, T- 
dependences of usp and toaf are plotted for zJ± = 0.7, 
when the staggered magnetic field is varied. With con- 
siderably small values of Hg, usp becomes completely 
suppressed, and at the same time toaf is drastically in- 
creased. The inflection point in the T-dependence of 
TO,AF J representing the characteristic temperature for the 
emergence of AF, is shifted about AT — 0.2 to higher-T, 
which is ten times the energy of the applied staggered 
magnetic field of Hg = 0.02. Thus we have demonstrated 
that there is actually a "colossal" response in this system. 

Now we compare our results to real systems. In doped 
CuGeOa , in some range of T above the impurity induced 
uniform AF phase, the static SP lattice distortions show 
a re-entrant behavior observed by elastic neutron scat- 
tering [18] as mentioned previously, similar to what is 
seen in our calculations (Fig. 4). However in experiments 
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FIG. 5. Temperature dependence of usp and mAP for 
zj±=0.7 under small staggered magnetic field, Hs- 



the correlation length remains finite [18], namely, the SP 
state remains to be SRO, while in our calculations it at- 
tains LRO. This discrepancy may be due to the treatment 
of the impurities in our calculations, such as to ignore ef- 
fects of the random impurity positions between neighbor- 
ing chains and/or to assume that the impurity modifies 
only Ki neglecting the change in the intrachain coupling 
J. Such effects would make the system more disturbed 
resulting in frustration for the SP state in the interchain 
directions as pointed out in Ref. [24], due to sohtons of 
SP order parameters changing the phase at the impuri- 
ties [5,6,9]. This actually reduces the correlation length of 
SP order parameter, but the SP lattice distortions would 
remain static as in the experiments even if such effects 
were considered. Therefore, we believe that main fea- 
tures of our results are relevant to doped CuGeOs. 

We note that another theoretical study in Ref. [7] also 
showed apparently similar re-entrant behavior in a nar- 
row parameter range. However, because of their weak- 
coupling continuum treatment with white noise-type dis- 
order, all of their states are with spatially uniform SP and 
AF order parameters. This is very different from the ex- 
periments [2,3,18], as well as from our SP and co-existent 
states, where we believe that the spatially inhomogeniety 
is crucial for the emergence of the huge response. 

Our observation that very weak perturbation to the 
system can drive such "colossal" response should be rel- 
evant to the actual compounds. Although the staggered 
magnetic field that we have applied is a fictitious field, 
perturbations that can change the relative energies of SP 
and AF states might also produce such drastic changes. 
For example, the application of pressure, p, is one can- 
didate, through the change in J_l/J or Ki resulting in a 



large change of the transition temperatures as we have 
seen in the phase diagram, Fig. 2. In doped CuGeOs, the 
hydrostatic pressure drives the system toward the stabil- 
ity of SP state, contrary to the usual expectations [19], 
due to the enhancement of the next-nearest-neighbor ex- 
change along the chains [2,25]. The enhanced dTc/dp in 
doped samples compared with the pure sample may be 
related to the sensitivities of our impurity induced states. 
Furthermore, an estimate of uniaxial dT^/dp^s in terms 
of the Ehrenfest relation using the measured specific heat 
and thermal expansion [26] also gives a large value. The 
uniform magnetic field might lead to the giant response 
as well. Up to now, the experiments are mainly done for 
small- a; samples which found no such large change while 
the system converts into an incommensurate phase [2,27]. 
Measurements in the multicritical regime are desired. 

In conclusion, we have demonstrated that a model of 
disordered SP system competing with the AF order shows 
a re-entrant behavior and huge responses to small per- 
turbations such as the staggered magnetic field or the 
change in the interchain coupling J±. This is analogous 
to the CMR effect in manganites, where the multicriti- 
cal phenomenon influenced by the disorder leads to the 
"colossal" response. 

We acknowledge Y. Tokura and Y. Tomioka for useful 
discussions and providing us data prior to publication, 
and M. Saito for critical readings of the manuscript. 
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